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; Abstract 

The resonant activation effect (RA) has been weil stud- 
ied in different ways during the last two decades, ft 
consists in the presence of a minimum in the mean time 
spent by a Brownian particle to exit from a potential 
well in the presence of a fluctuating external force, as a 
function of the mean frequency (or the correlation time) 
of the latter. This work studies the role played by the 
asymmetry of a pieccwise linear potential in the RA ef- 
fect, and, in general, the behavior of the mean first pas- 
sage time and the mean velocity of the particle crossing 
through the potential barrier. A strong dependence on 
the asymmetry of the potential has been found which 
can be put in relationship with the current in the ratchet 
whose the potential here used is an elementary module. 
In this case a current reversal as a function of the fre- 
quency of the switching potential occurs. Comparison 
of the calculations with the Doering-Gadua model have 
been performed, as well as comparison with smooth sym- 
metrical potentials, by checking for the robustness of the 
resonant correlation time. The calculations have been 
done by solving numerically the Langevin equation in 
the presence of an uncorrclatcd Gaussian noise. The 
resonant mean first passage times show an unexpected 
behavior as a function of the thermal noise intensity. 
The related curves present for the different symmetries 
an unexpected inversion of their relative behavior be- 
yond a certain threshold value of the noise. This means 
that the current reversal can only occur for weak noise 
intensities, lower than that threshold value. 

Pacs: 05.40.-a, 05.45.-a 

In the recent past years various theoretical works have 
been produced around the concept of Resonant Activa- 
tion (RA) [J i i B | Ej, which consists in the 
presence of a minimum of the mean escape time from 
a potential well of a Brownian particle when the sys- 
tem is subjected to a randomly switching force, as a 
function of the mean switching quantity. The RA effect 
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has been also detected experimentally [El 02 [TO] and it is 
in principle involved in a wide branches of science from 
physics to biology. The occurrence of the RA together 
with other stochastic effects such as noise enhanced sta- 
bility (NES) PU EH and stochastic resonance (SR) 
(T3] have been also investigated 02 [14] . The article by 
Doering & Gadua [T] is considered as one of the most in- 
troductory work to the resonant activation phenomenon. 
They introduced a switching piecewise linear potential in 
a range [0, L] with fixed minima in x = and x = L, 
and fluctuating amplitude of the maximum. They re- 
port the mean escape time for a Brownian particle in 
the case of fluctuations of the maximum of the potential 
between Vq and (flat potential), and also between Vq 
and —Vq (well instead of barrier), showing the presence 
of the RA effect. A slightly different choice was made by 
Bier & Astumiam [2] who used a potential fluctuating 
between Vq — a and Vq + a with a < Vq, maintaining so 
the presence of the barrier in all the dynamics. Both the 
methods show the resonant activation effect, and have 
been analytically evaluated in some approximation [Sj. 
Both the choices have in common that the potential is 
symmetrical in shape and maintains the same value at 
the two extrema in all the dynamics (V(0) = V(L)). 

Aim of this work is to focus on the role played by the 
asymmetry of the potential on the resonant activation 
effect using a simple piecewise linear potential . 

Many papers with both experimental and analytical 
investigation concerning the role played by the asym- 
metry of the potential in stochastic effects have been 
published during the last years. However, the investi- 
gations have been mainly devoted to the effects on the 
Stochastic resonance phenomenon, that is the presence 
of a noise induced regular oscillations in a system, which 
is revealed by means of a maximum in the signal to noise 
ratio of the output QH [HI TH UHl EES] , while the relation 
between the RA and the shape of the potential has been 
previously performed using a single slope linear potential 

US- 

Comparison with the Bier-Astumian model have been 
performed, as well as comparison with smooth symmet- 
rical polynomial potentials. The results obtained have 
been extended to the most elementary ratchet potential, 
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giving explanation of the current reversal there found as 
function of the correlation time of the external force. 

The fluctuating potential V±(x, t) is here given as the 
sum of a static potential V(x) (with, again, V(0) = 
V(L)) plus an additional time-dependent U(x,t) giv- 
ing the two configuration 'up' and 'down' between which 
V±(x, t) takes its values. The additional potential U(x, t) 
has not to be necessarily a stochastic process to give rise 
to RA [U[S]. It can be a smooth, continuous potential 
like a cosine or, instead, a stochastic potential related 
to a dichotomous force exponentially correlated in time. 
This last form is widely used in literature and we use it 
in this work. The related Langevin equation is: 



= -V'(x)+r)(t)+m 



(1) 



where is the Gaussian white noise, with zero mean 
and correlation function (£(*)£(*')) = 2D6(t-t'). The in- 
tensity D is related to thermal bath and damping coeffi- 
cient 7 (here 7 = 1) by means of the relation D = ^yksT. 
The random force rj(t) represents a dichotomous stochas- 
tic process, the random telegraph noise (RTN), taking 
the two values {—a, a} with an exponential correlation 
function (ri(t)r)(t')) = (Q/T)e~' t_ * '/ r , where the inten- 
sity Q = cl 2 t and t is the correlation time of the process. 
The potential V±(x, t) is then defined as: 
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Figure 1: Piecewise linear potential used in the calcu- 
lation. The position xq = represents the starting 
point of the simulations and there is also put a reflect- 
ing boundary. With respect to Bier-Doering choice (left) 
the right extremum of the potential is not fixed (center) 
and the whole potential flips randomly between the two 
shapes with a correlation time r. The right 

draw shows an example of asymmetric piecewise linear 
static potential. The parameter k represents the posi- 
tion of the maximum x m with respect to the position of 
the maximum in the symmetric case x s (here x s = 0.5, 
x rn = 0.25, k = -0.25) 



V±(x,t) = V(x) + U(x,t) =V(x) -xr)(t). (2) 



with, explicitly, 

V(x) = 



0<x r 



L-x 



(3) 



Here L = 1, h = 2, x m = L/2 + k, and k represents 
the asymmetry parameter, defined as the distance of the 
position of the maximum of the potential x m from the 
position of the symmetrical maximum x s . In Fig. [T] 



(right) we see an example of the static potential with 
the asymmetry parameter k = —0.25. 

The difference between the choice of the fluctuating 
potential here used (Eqs. [2] and [3]) and that one by Bier- 
Astumian is visible in Fig. [H where the potentials are 
drawn in the two cases. Here, with respect to the 'mean' 
static potential, the additional fluctuations —xrj(t) give 
only two values ('up' and 'down' in fact) in the flipping, 
being the force n(t) uniform overall the arrange of the 
potential. In the Doering & Gadua model (as in the 
Bier-Astumian) two values of the force for each potential 
slope have to be considered to hold the minimum on the 
right at the same level (V(L) = const). 

The potential V±(x,t) here defined can be considered 
as a base modulus of the piecewise linear ratchet sub- 
jected to RTN widely used in literature [2"T | \T2 \ |2"5]. 
The equation ([1]) has been solved numerically by using 

RA - Asymmetric PWL Potentials 
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Figure 2: MFPT showing the resonant activation effect 
for three values of the asymmetry parameter k of the 
piecewise linear potential (A; = —0.25, k = 0, such as 
symmetric potential, k = 0.25). The white noise inten- 
sity is here D = 0.18. The intensity of the dichotomous 
force is a = 1.2. The bottom/left inset show the mean 
velocity of the Brownian particle for the same asymme- 
tries, as a function of the correlation time r. 



dt = 10 -3 and the averages have been performed over a 
sampling of N = 20, 000 realizations. In the i-th realiza- 
tion the particle is put in the starting position xq = 
and the time ti to cross the position x = L is computed. 
A reflecting boundary is put in the left extremum of the 
potential while an absorbing boundary is present at the 
right extremum. The ensemble average of the ti gives the 
Mean First Passage Time (MFPT), which presents, for 
all the cases here studied, the evidence of the RA effect, 
i.e. a well drawn minimum as a function of the correla- 
tion time t. In fact, as well as the symmetric case, the 
MFPT obtained with the asymmetrical potentials show 
a resonant effect which is drawn in Fig. [2l We notice 
that for the three values of the asymmetry parameter fc, 
we find quite the same value of the resonant correlation 
time tr ~ 10, but different values of the corresponding 
resonant MFPTs (Trs), which decrease by increasing the 
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asymmetry parameter k. 

We note that the resonant region shows an inversion 
in the behavior of the MFPT curves for the three po- 
tentials to both the low and high correlation times with 
respect to the intermediate one. In fact for r lower than 
T c L ~ 10 _1 the curves show a MFPT higher for positive 
asymmetry (fc = 0.25) and lower for negative asymme- 
try (k = —0.25) and the same qualitative behavior is 
visible in the long correlation time region r higher than 
T c R ~ 10 3 - In the intermediate region r £ [ T c L -, T C R ], 
where wc also find the resonant values, the situation is 
inverted: the highest Tr value corresponds to the neg- 
ative asymmetry parameter and the lowest Tr to the 
positive one. 

On the other hand, calculation performed with the 
Bier-Astumian and Doering & Gadua model, that is us- 
ing fixed extrema of the same asymmetric piecewise lin- 
ear potentials, which fluctuates between the same highs 
(a = 0.6) in all the asymmetries, give strongly differ- 
ent curves, as visible in Fig. [31 where we can see even 
a strong displacement of the resonant correlation time 
by changing the asymmetry, but no crosses are present 
between the MFPT curves. The inversion of the MF- 

RA - PWL Pot. with fixed extrema 
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Figure 3: MFPT showing the resonant activation effect 
for different asymmetrical piecewise linear potential with 
fixed extrema and equal potential excursion in the three 
cases. Differently from the results in Fig. [5] we don't 
find any cross between the three curves and also the 
resonant correlation time tr changes for the different 
asymmetries. The potential increase is here a = 0.6, as 
the one of the symmetric case in Fig. [5] 



but, instead, because of the presence of the asymmetry 
together with the uniformity in space of the fluctuating 
external force rj(t) added to the system. 

The main relevant feature in adding a uniform force 
in the range of the constant potential, lies in the fact 
that in this case the barrier high of the fluctuating po- 
tential takes different values for different positions of the 
maximum, i.e. as a function of the asymmetry param- 
eter k. In fact the resonant MFPT values Trs depend 
mainly by the lower value taken by the potential (V-), 
being proportional to (l/V /2 )e l/ -/ 15 [Jj], and this value 
becomes lower and lower, by increasing the value of ks. 
This means that at a first sight we can expect that the 
Tr values take a lower value for the positive asymmetry 
than for the negative ones. However, as we can see below 
in the text, this expectation holds only up to a certain 
threshold value of noise intensity (Dt) and the inverse 
behavior occurs for higher values (D > Dt)- 

The model here investigated presents interesting fea- 
tures in the MFPT: first of all it has a value of the reso- 
nant mean period tr not too strongly dependent on the 
asymmetry parameter fc; then, it presents two period in- 
tervals, close to tc l ~ 10" 1 , and close to tc r ~ 10 3 
having approximatively the same MFPT for all the k- 
parameters. 
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Figure 4: Resonant activation evolution for various noise 
intensities. We observe the disappearing of the cross- 
ings in the MFPT curves and the shift of the minima 
by increasing the noise intensity D. These behaviors are 
shown in details Fig. [5l where the resonant frequency 
and resonant times are plotted as a function of D. 



PTs curves behavior, and consequently the presence of 
the two crosses at approximatively tc l and tc r , is so 
uniquely present in MFPTs calculated for asymmetri- 
cal potentials using uniform fluctuating force overall the 
range [0,L], and it does not appear neither in the sym- 
metrical ones with different shapes (See Fig. [7J, nor in 
the asymmetrical ones with fixed extrema and fluctu- 
ating barriers (Fig. [3]). In other words, the comparison 
between the results plotted in Figs. [21 [3] and [7] put in evi- 
dence that the cross features of the MFPT curves occurs 
not merely because of the asymmetry of the potentials, 



However, the crossing features of the MFPT as a func- 
tion of the mean driving frequency doesn't occur for any 
value of the thermal noise intensity D. A set of cal- 
culation to check this kind of robustness has been per- 
formed and the related results are shown in Fig. O where 
the RA is plotted for different values of the noise. We 
can see that by increasing the thermal noise intensity 
D, the crosses between the curves are maintained up 
to a threshold value that we can call Dt- For higher 
noise intensities no crosses appear in the curves. Fur- 
ther, a shift of the resonant mean switching time as a 
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function of noise is visible and a lowering of the related 
Tr, which demonstrate that when the noise is increased, 
an higher frequency switching is necessary to reach the 
resonance, and this resonance occurs at a lower mean 
escape time. The increase of the noise intensity has in 
this sense the effect to speed up all the escape features 
from the well. The results shown in Fig.[4j can be better 
observed in Fig. [5J where the mean resonant frequencies 
7r = 1/(2tr) and the mean resonant escape times Trs 
have been plotted as a function of the noise intensity for 
the three symmetry values. We can see there that go- 
ing beyond the threshold noise value Dt, the three Tr 
curves invert their relative position. This noise thresh- 
old corresponds to the presence (for D < Dt) or the 
absence (for D > Dt) of the two crossings of the MFPT 
curves visible in Fig. [2] and Fig. [U However we notice 
that Dt is not unique for all the asymmetries, being the 
crossing values slightly different for each couple of the 
three curves. We can see that the resonant frequency 
(inset of Fig. [5J has a slightly different dependence on 
the thermal noise intensity D for the different asymme- 
tries. In the range of D investigated, the three curves 
can be easily approximated by a straight line, even if the 
real dependence is in general more complicated (see [5]); 
the slope of this line is higher for the negative asymmetry 
than for the positive one. 

Resonant frequency and Mean first passage time vs D 
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Figure 5: Resonant mean passage times Trs as a function 
of the noise intensity D for the three asymmetry values 
investigated. In the inset the corresponding resonant fre- 
quencies 7rs. The crossing value Dt ~ 0.27 represents 
the threshold of thermal noise discriminating if the two 
crossings in MFPT of Fig. [5] are present (D < Dt) or 
absent (D > Dt). 



The presence of a resonant behavior, as well as the 
cross value at tq l , is also found in the plot of the mean 
velocity of the Brownian particle. Left inset of FigfS] 
shows, in fact, this measure as a function of the cor- 
relation time of the fluctuating dichotomous force, cal- 
culated as v = N^ 1 YliLi L/ti. For all the asymmetry 
parameters, we see the presence, before the saturating 
behavior, of a weak maximum which corresponds to the 
resonant correlation time tr. We can also see that for 



low values of the correlation times (r < Tc L ) the mean 
velocity is higher for negative asymmetry and lower for 
positive ones, while for (t > tc l ) is the inverse. This 
feature gives rise to a reversal current in the ratchet, as 
predicted in other works (2H [23 US] and whose occur- 
rence has been also demonstrated experimentally [2"7] . 
In fact the difference between the mean velocities of the 
positive asymmetry and the negative one change sign 
at the tc l value. In an asymmetrical ratchet this dif- 
ference represents a net velocity flux, provided that the 
absence of any reflecting boundary in that case gives rise 
to changes in the values of the velocity. Both the pres- 
ences of a maximum for r ks tr and the cross at r ps tc l 
are in total agreement with the behavior of the MFPT. 
This agreement fails, instead, for values of the correlation 
times higher than tr. While the MFPT curves increase 
in a different way and joint together at the second cross, 
the velocities decrease only a few, reaching a saturation 
value. This is because for high values of the correlation 
times, the particle tends to cross the potential barrier 
when it is in its lower high, so acquiring a relatively high 
speed because of the low travelling time. When the po- 
tential is in the high level, the particle takes a longer 
time to cross and so the contribution to the mean ve- 
locity becomes very low and relatively negligible. This 
means that, for high correlation times, v maintains a 
relatively high value which doesn't change so strongly as 
the MFPT. 

RA - Elementary Ratchet 
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Figure 6: MFPT showing the resonant activation effect 
for the simplest ratchet potential, composed by only two 
the single barrier. The MFPT is here the mean time 
taken by the Brownian particle starting at x = to 
reach the position x = 1 or x = — 1. The mean velocity 
of the particle is plotted in the upper/left inset. The 
mean velocity has again a maximum at the same reso- 
nant value tr. In the very low correlation time region, 
the system has a weak negative velocity (right/top in- 
set), as a consequence of the different behavior and the 
cross of the mean velocities shown in the inset of Fig. [5] 



The results found above for the single barrier poten- 
tials, mirrors, of course, to the ratchet potential having 
the same asymmetric profile as elementary module. In 
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this respect a set of calculations has been performed with 
the aim to join together the results of the single barrier 
described above with the simplest ratchet case, such as a 
ratchet with two barriers only. Fig. [H] shows the results 
in such a case and the bottom/right inset shows the cor- 
responding elementary ratchet. The system consists of 
two asymmetric barriers without the presence of any re- 
flecting boundary. The MFPT presents again a resonant 
correlation value tr which is the same for the single bar- 
rier case, as we can expect. In this system the MFPT is 
the mean time spent by the Brownian particle starting 
at x = to reach the position x = 1 or x = —1, in- 
differently. The particle, of course will follow the easiest 
path, and the MFPT represents the minimum time of the 
two single barrier case seen above. This also means that 
the curve is lowered and the RA effect less pronounced. 
The mean velocity, plotted in the upper-left inset of the 
Fig. [6j shows again a maximum at the same resonant 
value tr. For very low correlation time the mean ve- 
locity has a weak negative velocity (right /top inset in 
FigE]). This means that a current reversal appears at a 
certain correlation time r rev . This features follows from 
the different behavior of the mean velocity in the two 
specular asymmetric single barrier potentials seen above 
(inset of FigJ^J), where the presence of the cross value 
t<7 £ indicates a current reversal as a function of r. The 
difference in value between r rev and tq l , as well as the 
difference in the absolute value of the mean velocity of 
the Brownian particle, have to be imputed to the pres- 
ence of the reflecting boundary in the single barrier case 
which change the traveling times of the particle and, so, 
the related mean velocities. 

RA - Various potential shapes 
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Figure 7: MFPT showing the RA effect for different 
shape of symmetrical potential. We note than for r < tr 
the logarithmic distances between two MFPTs is approx- 
imatively a constant. This evidence suggests that in that 
region an exponential from factor should be inserted in 
the analytical expression of the mean first passage time 
to take into account the shape of different potentials. 
The parameters arc the same than in FiglU 

As a last remark concerning the relationship between 
the resonant activation effect and the shape of the poten- 



tial, some calculations have been performed using sym- 
metrical smooth potentials. The static potentials used 
have the form: 

V 2N (x) = h2^(l-^) (4) 

where in our calculation h = 2. The values used are: 
N = 1,2,3, such as parabolic, quartic and 6 th power 
potentials. As we can see in Fig. [7J the resonant mean 
time is quite the same for all the cases, again confirm- 
ing that tr is a robust value in the model investigated. 
Another remarkable and well visible feature is that the 
four curves of the MFPT differ each other of a constant 
quantity, at least for r < tr. This means that, in that 
region, their logarithmic distance is constant and so an 
exponential form factor has to be taken into account in 
order to estimate the MFPT for each potential shape. 

Summarizing the results, the shapes of the potential 
(both symmetrical and asymmetrical ones) play a very 
important role in the evaluation of the RA effect and 
MFPT behaviors. With a spatially uniform random tele- 
graph force, the resonant correlation time tr appears to 
be a robust value independently on that shape, while this 
latter acts always in a strong way by modifying the res- 
onant values of the mean first passage times Trs. In the 
context of uniform forces, the asymmetry of the potential 
is then responsible for the crosses of the MFPT curves in 
a certain range of low thermal noise intensities, giving an 
explanation for the appearance of the current reversal as 
a function of the correlation time of the fluctuating force 
in ratchet potentials. These crosses, and, consequently 
the current reversal in ratchet, are only present at weak 
noise intensity, as indicated by the presence of an upper 
noise intensity threshold Dt- 
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